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Using an advanced tight-binding approach, we estimate the anisotropy of the tunnel transmission 
associated with the rotation of the 5/2 spin of a single Mn atom forming an acceptor state in 
GaAs and located near an AlGaAs tunnel barrier. Significant anisotropics in both in-plane and 
out-of-plane geometries are found, resulting from the combination of the large spin-orbit coupling 
associated with the p-d exchange interaction, cubic anisotropy of heavy-hole dispersion and the low 
C2v symmetry of the chemical bonds. 
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In the context of semiconductor-based spintronics 
GaMnAs-based tunnel diodes are of utmost interest. In 
particular, it was recently demonstrated that the tunnel- 
ing current in semiconductor heterostructures integrat- 
ing the p-type ferromagnetic semiconductor GaMnAs is 
strongly affected by the direction of its magnetization, 
which can be arbitrarily changed by an external magnetic 
field 1-4 . In this regard, two main phenomena must be dis- 
tinguished: (i) the tunneling magnetoresistance (TMR) 
which is the dependence of the tunnel current with the re- 
spective parallel or antiparallel magnetic configurations 
of two ferromagnetic layers separated by a tunnel bar- 
rier giving rise to spin-valve effects and (ii) the tunnel- 
ing anisotropic magnetoresistance (TAMR) which is the 
change of the tunnel current as a function of the magneti- 
zation direction of a (single) magnetic layer with respect 
to a reference coordinate system. So far, experimental 
results were obtained in the metallic regime, but it 
was suggested that the origin of TAMR could be still im- 
puted to the anisotropic shape of the hole state bound to 
the isolated Mn impurity in a GaAs host . Various the- 
oretical approaches indicate that the bound hole state 
is highly anisotropic and more extended in directions 
perpendicular to the Mn spin orientation. As a result, 
for in-plane magnetization, the exponential tail of hole 
bound state is more extended and tunneling probability 
higher. On the one hand, the validity of this explanation 
might be a matter of discussion for Mn densities close 
to the Mott limit (few percent Mn, considering a Bohr 
radius of 1 nm). On the other hand, even this qualitative 
picture needs to be checked by calculations including the 



effect of the barrier on the impurity state itself. Indeed, 
it is known from STM studies that hole wavefunction is 
strongly affected by its structural environment 8-10 . In 
this Letter we consider a model system close to experi- 
mentally achievable single impurity tunnel diode and es- 
timate the TAMR effect in the framework of an advanced 
sp 3 d 5 s* tight-binding model 11 including exchange inter- 
actions in the effective semi-classical approach 6,7 . In ad- 
dition to a strong dependency of tunnel probability on the 
angle between Mn spin orientation and the [001] growth 
direction, our calculations reveal an unexpected in-plane 
anisotropy of comparable magnitude, with the [110] and 
[110] eigenaxis characteristic of the C^v symmetry 12 ' 13 . 

We model TAMR using the following scheme: we 
consider a single Mn impurity in a GaAs layer near 
a thin AlGaAs barrier (see Fig. 1) and calculate the 
neutral acceptor bound state. Tunneling transmission 
from this hole bound state through the AlGaAs barrier 
is estimated as an integral of the wavefunction tail 
in the GaAs "collector" layer on the right side of the 
barrier (z > 0): (^(r)) J z>Q \^h(f)\ 2 d 3 r . We do not 
consider here the "injector" problem, i.e. the way a hole 
is injected onto the impurity bound state from the left 
of the structure. This step can be ensured by p-doping 
the left electrode with shallow acceptors like Be, as was 
done in the context of STM studies . 

For the physical validity of our calculations, all what 
matters is that the current be limited by the tunnel bar- 
rier in Fig. 1. The hole density of states (DOS) at the 
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FIG. 1. Schematic representation of a Mn impurity in the 
GaAs layer near a thin Alo.3Gao.7As barrier. Wavefunction of 
the hole bound to the acceptor penetrates this barrier. Tun- 
neling probability is proportional to the integral of the hole 
wavefunction over the region to the right of the barrier. 
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right collector like originally investigated by Bardeen 14 
is not explicitely included, as it just enters a prefactor 
that does not depend on Mn spin orientation. Still, this 
simple approach is only a zeroth order estimate of the 
tunnel probability as it does not take into account the 
spin structure of hole wave- functions. One may convince 
himself that this should be valid as soon as the kinetic 
energy at which the hole is injected into the collector 
is high enough (a few 10 meV's). Indeed, in this case 
the admixture of heavy and light characters is such that 
the selectivity of transmission into the heavy- and ligh- 
hole bands should essentially vanish. A rough estimation 
shows that this corresponds to actual experimental situa- 
tion. Note that a more complete treatment would require 
full description of a specific collector geometry. 

To calculate wavefunctions we use the model and pa- 
rameterization of Ref. 9 and include p-d exchange inter- 
action in the mean-field approach previously introduced 
by Tang and Flatte . This well established method, valid 
in the ferromagnetic regime, actually corresponds to an 
Ising coupling between the hole and Mn 3d 5 electron 
spins. It amounts to treating the Mn 5/2 spin as a clas- 
sical magnetic momentum . In the paramagnetic regime 
a complete quantum treatment based on Heisenberg ex- 
change coupling would be definitely required, but it 
is much more difficult to implement self-consistently in 
the tight-binding formalism. It is also worth mentioning 
that for a given (classical) Mn spin orientation, the mag- 
netic acceptor ground state is non-degenerate, in contrast 
with the four-fold degeneracy of a non-magnetic neutral 
acceptor state. 

In practice, we use a supercell formalism, with a 8624 
atom supercell of 7 MLx7 MLx22 ML, which is approx- 
imately 4 nmx4 nm in lateral directions and 12 nm in 
growth direction. Technically, we destroy the effects of 
periodicity in the [001] growth direction by adding an 
AlAs barrier to the left of the GaMnAs layer, that fully 
decouples adjacent cells. Periodic boundary conditions 
in lateral direction produce non-negligible effects on the 
bound state energy and on the tails of the wavefunction, 
that must be carefully characterized. However, these ef- 
fects may hardly be considered as an artefact of the model 
as long as distance between Mn atoms corresponds to less 



FIG. 2. (Color online) Cross-section of the neutral acceptor 
wavefunction in the (110) plane containing the Mn dopant. 
Wavefunction amplitude is color-coded according to the log 
scale. Atomic positions are shown with circles : red and 
blue for As and Ga, green for Mn, and white for the virtual 
Alo.3Gao.7 cation in the tunnel barrier. 



than one Mn atom per thousand in the GaMnAs alloy. 
In real structures Mn atoms would be distributed ran- 
domly in the alloy while in the calculations they form a 
periodic array on a square lattice. We have checked that 
calculation results are stable with respect to Mn spatial 
distribution. Another point worth to be mentioned is the 
use of sp 3 d 5 s* model, compared to the sp 3 model used 
in previous papers ' . While sp 3 is a qualitatively ac- 
ceptable basis to discuss valence band properties in III- 
V compounds, it is well-known that it fails to repro- 
duce accurately the Luttinger parameters that describe 
valence band dispersion. sp 3 d 5 s* is actually the minimal 
tight-binding model which allows fitting all the param- 
eters of valence band over a large energy range, and in 
particular, the cubic anisotropy (warping) of the heavy- 
hole dispersion. This is essential for an accurate account 
of the kinetic energy part of the impurity hamiltonian. 
In order to image the wavefunctions, the local density of 
states at atomic sites is approximated by gaussian func- 
tions with a radius depending on orbital type (~1.5-2A). 

The cross-sections of the wavefunction of a hole bound 
to the Mn acceptor are shown in Fig. 2 for two direc- 
tions of the Mn spin, in the (llO)-plane containing the 
impurity. The wavefunction amplitude is color-coded ac- 
cording to the log scale below the panels, and atomic 
positions are indicated by colored dots (see figure cap- 
tion). Left panel displays the situation for a Mn spin 
along the [001] direction, and right panel for a Mn spin 
in [110] direction. While it is not obvious from these fig- 
ures how tunneling probability can be affected, it may 
be concluded that the hole wavefunction changes signifi- 
cantly in size as already emphasized in Ref. . The effect 
of barrier proximity is clearly visible from the left /right 
asymmetry of wavefunction tails, and the role of lateral 
periodicity is evident from the non vanishing tails in the 
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FIG. 3. Dependence of wavefunction integral on right side 
of barrier on Mn spin orientation, a) dependence on zenithal 
angle for different values of the azimuthal angle ip with 
respect to [100], b) azimuthal plot of in-plane anisotropy. 



vertical direction in Fig. 2 (supercell in the calculations is 
few ML larger in lateral direction than shown in figure). 
Neutral acceptor binding energy decreases when the im- 
purity gets closer to the barrier . For Mn located at 1 
ML from the barrier, this reduction amounts to 20 meV. 
The eigenenergies also change with Mn spin orientation: 
this is associated with the magnetization-dependence of 
the wavefunction tails, and the related change in tunnel 
coupling of acceptors in adjacent cells . The amplitude 
of this change for Mn located at 1 ML from the barrier 
as a function of angle between Mn spin and growth di- 
rection is almost 10 meV. Under in-plane rotation of the 
spin direction this modulation of eigenenergies is reduced 
but still amounts to about 4 meV. The calculated tun- 
neling probability is extremely sensitive to the impurity 
distance from the barrier, decreasing exponentially with 
characteristic length ~3.7 A. However, for all distances 
it shows similar angle dependencies. In the following, we 
focus on the 1 ML distance corresponding to Fig. 2. 

Calculation results are illustrated in Fig. 3. Mn spin 
direction is defined with the zenithal angle with re- 
spect to [001] growth direction and in-plane azimuthal 
angle ip with respect to [100] axis. The dependency on 6 
follows the qualitative explanation given in the introduc- 
tion: the wavefunction of the hole localized at acceptor 
has a somewhat oblate shape and follows Mn spin direc- 



tion. If the spin is oriented along the growth direction, 
the barrier lies in the direction where the exponential de- 
cay is stronger, giving a lower tunneling probability by 
more than 20%, which was defined in a previous letter 
as a negative TAMR . Wavefunction deformation due 
to the barrier proximity and large cubic anisotropy due 
to valence band warping (see Fig. 2) do not change this 
qualitative argumentation significantly. In-plane varia- 
tion (inset in Fig. 3) is more intriguing. It has similar 

amplitude as the zenithal variation, and shows the char- 
acteristic Civ symmetry, with pronounced differences be- 
tween [110] and [110] directions. We have checked that 
this result is not a calculation artifact due to lateral peri- 
odic boundary conditions by rotating the supercell orien- 
tation, which is equivalent to rotating the square lattice 
of Mn atoms: this rotation has no influence on the calcu- 
lated tunnel current anisotropy. This suggests that these 
results would be qualitatively unchanged for the ran- 
domly positioned Mn atoms. In-plane anisotropy is also 
nearly unaffected by the lateral supercell size, although 
the variation of bound-state eigenenergy with spin ori- 
entation increases with decreasing size (or, equivalently, 
with increasing Mn concentration). The calculated in- 
plane anisotropy of about 14% actually results from a 
combination of three factors reducing the symmetry of 
bound state wavefunction : i) a native anisotropy of the 
bulk wavefunction, which is clearly evidenced by remov- 
ing the barriers in the supercell while calculating the in- 
tegral over the same part of the wavefunction ; ii) the 
left /right asymmetry of the wavefunction due to the in- 
terface proximity ; and iii) the interface Civ symmetry 
itself ' . These contributions interfere and can hardly 
be separated in the final result. 

In summary, we have estimated the anisotropy of the 
tunnel transmission from the fundamental A hole state 
attached to a single Mn dopant in a GaAs host matrix, 
coupled to a reservoir through an AlGaAs tunnel barrier. 
The simulated structure displays pronounced anisotropy 
of the tunnel transmission as a function of Mn classical 
spin direction. Significant in-plane anisotropy of the tun- 
nel current is obtained. Yet, a full treatment of Heisen- 
berg (instead of Ising) spin coupling and direct account 
of applied magnetic field would be highly desirable in or- 
der to fully assert the potential of single impurity TAMR 
effect. 

Acknowledgment: This work was supported by "Trian- 
gle de la Physique", by CNRS-RAS joint laboratory IL- 
NACS and by the European Networks for Initial Training 
(ITN): SemiSpinNet. 



1 C. Gould, C. Riister, T. Jungwirth, E. Girgis, G. M. 
Schott, R. Giraud, K. Brunner, G. Schmidt, and L. W. 
Molenkamp, Phys. Rev. Lett. 93, 117203 (2004). 

2 C. Riister, C. Gould, T. Jungwirth, J. Sinova, G. M. 



Schott, R. Giraud, K. Brunner, G. Schmidt, and L. W. 
Molenkamp, Phys. Rev. Lett. 94, 027203 (2005). 
3 K. Pappert, M. J. Schmidt, S. Hiimpfner, C. Riister, G. M. 
Schott, K. Brunner, C. Gould, G. Schmidt, and L. W. 



4 



Molenkamp, Phys. Rev. Lett. 97, 186402 (2006). 

4 E. Y. Tsymbal and I. Zutic, eds., Handbook of Spin Trans- 
port and Magnetism (Chapman and Hall/CRC, 2011) p. 
808. 

5 M. J. Schmidt, K. Pappert, C. Gould, G. Schmidt, R. Op- 
permann, and L. W. Molenkamp, Phys. Rev. B 76, 035204 
(2007). 

6 J.-M. Tang and M. E. Flatte, Phys. Rev. B 72, 161315 
(2005). 

7 T. O. Strandberg, C. M. Canali, and A. H. MacDonald, 
Phys. Rev. B 80, 024425 (2009). 

8 P. M. Koenraad and M. E. Flatte, Nat Mater 10, 91 (2011). 

9 J.-M. Jancu, J.-C. Girard, M. O. Nestoklon, A. Lemaitre, 
F. Glas, Z. Z. Wang, and P. Voisin, Phys. Rev. Lett. 101, 
196801 (2008). 



10 O. Krebs, E. Benjamin, and A. Lemaitre, Phys. Rev. B 
80, 165315 (2009). 

11 J.-M. Jancu, R. Scholz, F. Beltram, and F. Bassani, Phys. 
Rev. B 57, 6493 (1998). 

12 O. Krebs and P. Voisin, Phys. Rev. Lett. 77, 1829 (1996). 

13 E. L. Ivchenko, A. Y. Kaminski, and U. Rossler, Phys. 
Rev. B 54, 5852 (1996). 

14 J. Bardeen, Phys. Rev. Lett. 6, 57 (1961). 

15 A. K. Bhattacharjee and C. B. a la Guillaume, Solid State 
Commun. 113, 17 (2000). 

16 T. B. Boykin, G. Klimeck, R. C. Bowen, and R. Lake, 
Phys. Rev. B 56, 4102 (1997). 

17 G. Bastard, Phys. Rev. B 24, 4714 (1981). 

18 M. Elsen, H. Jaffres, R. Mattana, M. Tran, J.-M. George, 
A. Miard, and A. Lemaitre, Phys. Rev. Lett. 99, 127203 
(2007). 



